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Analytical relationships for the surface and curvature energies of oblate and prolate semi- 
spheroidal atomic clusters have been obtained. By modifying the cluster shape from a spheroid 
to a semi-spheroid (including the flat surface of the end cup) the most stable shape was changed 
from a sphere to a superdeformed prolate semi-spheroid. Potential energy surfaces vs. deformation 
and the number of atoms, A'^, illustrate this property independent of A''. 

PACS numbers: 36.40.Qv, 71.45.-d, 61.46.Bc 



The density functional theory [ij is successfully cm- 
ployed in the field of atomic cluster physics. Alterna- 
tiyely, with less computational effort, one can use as a 
first approximation some simple models (see the reviews 
and the references therein) replacing the many-body 
effects by an effective single-particle potential, since to a 
good approximation the delocalizcd conduction electrons 
of neutral small metallic clusters form a Fermi liquid like 
the atomic nucleus 

The hquid drop model (LDM) @ dominated for many 
decades the theory of nuclear fission, starting with the 
first explanation, given in 1939 by Lise Meitner and O. 
Frisch, of the induced fission process discovered by O. 
Hahn and F. Strassmann. One of the most cited paper 
in the field was published in the same year Q . We have 
used the LDM to develop the analytical superasymmetric 
fission model allowing to predict in 1980 heavy particle 
radioactivity 

In 1990 W. A. Saunders [i| adapted the LDM to the 
atomic cluster physics, and explained the increase of fis- 
sionability with decreasing size of the charged metal clus- 
ter which was observed in experiments. The following 
year J. P. Perdew et al. presented a LDM for a neutral 
metal cluster with z = 1,2,3,4 valence electrons. They 
mentioned in the abstract that the LDM "originally de- 
veloped for finite systems (nuclei), may actually be more 
appropiate for infinite ones (metals)." For the ground 
state properties of neutral clusters or the fission of doubly 
or multiply charged clusters, the LDM 0, [HI, [H, [H, 
expresses the smooth part of the total energy to which 
the shell corrections [3, 0, [13, [3 [3, [13, Ell, S may 
be added. The interplay between LDM, shape deforma- 
tions and Strutinsky shell corrections (including the case 
of fission) for free clusters was reviewed by C. Yannouleas 
et al 



An interesting application of the liquid drop surface 
tension for producing micro or nanoscale objects was re- 
cently reported We beheve that LDM completed 
with shell corrections [l3| could be very useful for the 



physics of the atomic clusters deposited on different types 
of surfaces. It may provide very rapidly a first solution 
which can be eventually refined within density functional 
theory. 

The nanostructured coating of surfaces by cluster de- 
position 24, 2^ is at present a rapidly growing field. 



25[ is at present a rapidly growing field. By 
analyzing some shapes of cluster deposited on a surface 
obtained by using scanning probe microscopy [H, \2^ . 
one can see that a semi-spheroid with the z axis of cylin- 
drical symmetry oriented perpendicularly on the surface 
plane may be a good approximation. We present ana- 
lytical results for the deformation-dependent surface and 
curvature energies, which has the advantage of shortest 
computer time and easiest interpretation. 

In all LDM studies published until now, both in nu- 
clear and atomic cluster physics, the most stable shape 
(minimum deformation energy) was a sphere. One would 
expect that for a semi-spheroid the most stable shape 
would be a semi-sphere or an oblate shape. For the first 
time, we found the surprising result that the superde- 
formed prolate semi-spheroid is in this case the equilib- 
rium shape. 

We investigate the stability of semi-spheroidal shapes 
by assuming, as a first approximation of one possibility 
which can be met in practice, a vanishing interaction en- 
ergy with the surface on which the cluster is deposited, 
so that the neutral atomic cluster may be considered to 
be free. Other types of shapes obtained from a spheroid 
by removing more or less than its half [2^ will be con- 
sidered in the future; the interaction with surface will be 
taken into account as well. Our choice is motivated by 
the fact that the corresponding shell model [2^ allows 
to obtain analytical formulae for the single-particle ener- 
gies. The remarkable result of this shell model is that for 
the first time the maximum degeneracy is associated to 
the superdeformed prolate shape. In this way both LDM 
and the shell structure suggest the enhanced stability of 
such a distorted shape. 

We are using the standard notation of (p, z) for the 
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axially symmetric dimensionlcss cylindrical coordinates. 
When the shape is a semi-spheroid the length scale is 
given by the radius of a sphere with the same volume, 
Rs = = 2i/3rsiVi/3^ in which N is the number 

of atoms, is the Wigner-Seitz radius (2.117 A for Na 



isl . 20l . 30j ) and p = p{z) is the surface equation given 
by 
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where a is the minor (major) semiaxis for prolate (oblate) 
semi-spheroid and c is the major (minor) semiaxis for 
prolate (oblate) semi-spheroid. Volume conservation 
leads to a^c = 1. 



It is convenient to choose the deformation 17| param- 
eter 6 defined by 
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so that 



a 
c 



(2) 



(3) 



The eccentricity is defined by the equation 
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FIG. 1: Binding energy per atom {—Eb/N) versus the num- 
ber of atoms A*' for Na clusters. Comparison between spherical 
and semi-spherical shapes. 
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The deformation energy with respect to semi-spherical 
shape of a neutral cluster is expressed as 

E-E''^ = {Es - Ef) + {Ec - Ef) 

= Ef[§p-l)+Ef(§p-l) (5) 

E-E^' = Ef {B-l,,,j - 1) + Ef {Bl^^, 1) (6) 
where for a semi-sphere one has 

Ef = SnRja = SnA^/^Rla = ^ijO (7) 



Ef = 2TTRs"fc = 2tt2^/^R, 
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with E'^ and -Ecuru the surface and curvature energy of a 
sphere with the same volume. When the atom is electri- 
cally charged, one has to add the Coulomb energy. 

Clusters and nuclei are "leptodermous" systems char- 
acterized by a constant density in the volume and a thin 
surface layer allowing to expand their binding energy in 
terms of powers of N^^"^ (see a detailed discussion in refs. 
0, IsH). Despite the fact that this expansion is a priori 
valid only for large enough systems (e.g. Na2654) it 
"predicts the energy per electron E/N accurately (within 
0.03 eV) even for N = 1". For a spherical Na cluster 
[isl . [2^ the binding energy, in eV, is given by 

En = -2.2b2N + 0.541A^2/3 ^ o.l54iV^/^ (9) 

where = —2.252N eV is proportional to the volume 
(assumed to be conserved), = 0.541A^^/'^ eV is pro- 
portional to the surface area and the surface tension a 



E"^ = AirRla = Airr^aN^^^ 



(10) 



hence Anria — 0.541 cV for Na clusters. The curvature 
energy E^^^^ = 0.154A^^/'^ eV is proportional to the in- 
tegrated curvature ant to the curvature tension 7c 



47ri?o7c = 47rr^7ciVi/3 



(11) 



where Airrgjc = 0.154 eV for Na clusters. 

The numerical coefficients in eq. ([9]) have been de- 
termined [13, by fitting the extended Thomas-Fermi 
local density approximation total energy [s^l for spheri- 
cal shapes. In the fig. 2 of ref. [31 the smooth line ex- 
pressed by eq. ^ is compared to the dots from ref. [3^ . 
The two sets of data coincide at magic numbers; shell ef- 
fects, explaining the deformation energy of non-spherical 
atomic clusters, may be added by using Strutinsky's [l5l | 
procedure. 

The liquid drop part (volume, surface, and curvature 
terms) of the binding energy of Na semi-spherical clusters 
will be in eV: 
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A cluster with a spherical shape is more tightly bound 
than a cluster with a semi-spherical shape, as shown in 
figure [TJ 

The area of a surface of revolution about z axis is 



y"dS' = 27ri?2 J p(^z)^i+(^^^ dz (13) 

The deformation-dependent surface energy for cylindrical 
symmetry Bsurf — S/Sq — S/AttRq is given by 
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dzp^l + p'2 = - 
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which for spheroidal (oblate or prolate) shapes becomes hence 



= a/(2c2) J^^ dzy/ci + z^ia^ - c^) (15) 
because for a spheroid 
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is 47ri?o- T^i and 7?.2 are the two principal radii of curva- 
ture at a local point on the surface. The shape-dependent 



The local curvature 

K = o.5(7^];^ + 7^^^) 

of a sphere is l/i?o and the integrated curvature 
K = I dSn 
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FIG. 2: Absolute values of the deformation-dependent liq- 
uid drop surface and total (surface plus curvature) energy 
versus the deformation parameter S for Nase clusters. Com- 
parison between spheroidal and semi-spheroidal shapes. The 
volume energy _E„ = —126.1 eV was not included in Eld and 
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part of the curvature energy for cylindrical symmetry is 
given by 
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dSK ^ / dzKp^l + p'2 (20) 
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The principal radii of curvature of a shape with axial 
symmetry 3^ are expressed as 

7^l = i?oPv/l + ^2 = -Ro- ^^jT— (21) 
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For the relative surface energy of a semi-sphere we ob- 
tain 

BsurfiS = 0) = (47ri?2/2 + 7ri?2)/(47ri?2) 

= (3/4)i?2/i?2 = 3/42/3 (23) 

and the corresponding curvature energy 

5c„™(0) = K/AttRo - (4^i?,,/2)/(4^i?o) 

= (l/2)(i?,/i?o) = 1/41/3 (24) 

We give p,z,a,c in units of Rs = 2^/^Ro, hence ac- 
cording to the volume conservation, a^cRg/2 — Rq so 
that a^c = 1 
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because for a semi-spheroid there is a plane circular sur- 
face with the area ttR^o^. The deformation dependent 
curvature energy of a semi-spheroid, B^^rv — K/{2nRs) 
is 
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dSK = -^ / dzK^p^ + {pp'f (27) 



The integrated curvature for the plane surface is zero, 
hence 



^curv — 2 Jo \ J 

2 ^ 2 Jo c4+z2(a2-c2) 
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and the second term of this equation can be simplified 
by taking into account that a^c = 1. 

When a > c (oblate semi-spheroid) and e'^ = a2/c2 — 1, 
a2 — c2 = e2c2, and 
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For a semi-spherical shape 

Asymptotically for a » c one has e ^ a/c and conse- 
quently 
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FIG. 3: Liquid drop deformation energy contour plot for 
Na clusters vs. deformation 5 and tiie number of atoms, N. 
Deformation energy relative to the semi-spherical shape. 

When c > a (prolate semi-spheroid), from the defini- 
tion of the eccentricity we have — = e^c^, hence 
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Similar equations for spheroidal shapes may be found in 
Ref. [s^l, where the notation 77 = a/c is used. 

Asymptotically for c ^ a one has e — > 1 and conse- 
quently 
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A comparison between the two LDM deformation ener- 
gies in figure [2] illustrates the essential difference between 
the spheroidal and semi-spheroidal shapes of a Na cluster 
with 56 atoms; in the former case the most stable config- 
uration is a sphere, but in the latter it is a superdeformcd 
prolate semi-sphcroid with S = 0.65 (c/a = 1.96). 



The contour plot of deformation energy vs. the de- 
formation (5 and the number of atoms N from figure [3] 
shows the general trend of stabilty of superdeformed pro- 
late shapes independent of the number of atoms in the 
cluster. 

Consequently by modifying the shape from a spheroid 
to a semi-spheroid the most stable shape was changed 
from a sphere to a superdeformed prolate semi-spheroid. 
In applications it will not be easy to measure the inter- 
action energy of the deposited cluster with the surface. 
The geometry may differ according to the strength of 
interaction between the cluster and the substrate. By 
observing in the experiment a prolate semi-spheroidal- 
like shape one can conclude that the interaction energy 
is very small. 
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